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Abstract — Wireless ad hoc networks are seldom characterized 
by one single performance metric, yet the current literature 
lacks a flexible framework to assist in characterizing the design 
tradeoffs in such networks. In this work, we address this problem 
by proposing a new modeling framework for routing in ad hoc 
networks, which used in conjunction with metaheuristic multi- 
objective search algorithms, will result in a better understanding 
of network behavior and performance when multiple criteria 
are relevant. Our approach is to take a holistic view of the 
network that captures the cross-interactions among interference 
management techniques implemented at various layers of the 
protocol stack. The resulting framework is a complex multiobjec- 
tive optimization problem that can be efficiently solved through 
existing multiobjective search techniques. In this contribution, we 
present the Pareto optimal sets for an example sensor network 
when delay, robustness and energy are considered. The aim of 
this paper is to present the framework and hence for conciseness 
purposes, the multiobjective optimization search is not developed 
herein. 



I. Introduction 

Wireless ad hoc or sensor networks often operate in difficult 
environments and require several performance criteria to be 
satisfied, related to timely, reliable, and secure information 
transfer To ensure information transfer across a network, 
one of the key elements is the selected routing protocol 
whose design poses significant challenges. In such networks, 
cooperation among all layers of the protocol stack should be 
enlisted to deal with channel impairments, and thus the design 
of a routing protocol should be viewed in the context of its 
interactions with other interference management techniques 
implemented at other layers of the protocol stack. 

To further add to the list of design challenges, it is seldom 
possible to "equally optimize" all desirable performance cri- 
teria, as some of them may be antagonistic in nature. From 
a myriad of possible operating points, which one is "more 
optimal"? Understanding the tradeoffs involved with respect 
to various performance metrics will not only lead to a better 
design, but also will allow for the selection of a set of possible 
operating points (characterized by various tradeoffs) to enable 
a graceful degradation of the network performance as the 
channel conditions worsen. While significant work has been 
done for routing in wireless ad hoc or sensor networks, no 
integrated design framework exists to address the many facets 
of the problem described above. 



There is a significant effort to characterize the theoretical 
performance of ad hoc wireless networks. Most of it is focused 
on their theoretical capacity, which has been assessed by 
several landmark papers under various assumptions |T|-f4l. 
However, none of these works directly supports a practical 
implementation of a routing algorithm, and they lack a general 
view of multiple objective tradeoffs - though some of them do 
consider the impact of the end-to-end delay on capacity. On 
the other hand, there is vast literature on designing routing 
protocols optimized for various specific criteria and specific 
network instances (e.g. |[5), |[6| and the references within). It 
is very hard to compare the quality of these solutions as no 
benchmarks for multiple criteria performance routing exist. 
Limited work exists on designing multiobjective (MO) routing 
1 7 1, and again the network scenarios used for optimization are 
very application specific. 

Finally, for networks operating in harsh environments, char- 
acterizing cross-layer interactions is essential as for instance 
in implementing interference mitigation at all layers of the 
protocol stack. Several seminal works have been written on 
various aspects of cross-layer design (e.g. (j3), (8)). 

Understanding the tradeoffs involved with various routing 
solutions will enable adaptive resource management across 
layers and nodes, leading to a more accurate "local to global 
performance mapping" for practical routing protocol design. 
Our main contributions in this work are two-fold: 

• Propose a general cross-layer framework network model, 
capable of capturing the impact and interaction of a 
wide range of interference and resource management 
techniques for various channel conditions; 

• Formulate a multiobjective routing optimization problem 
by defining appropriate evaluation functions for criteria 
such as: robustness of information transfer, end-to-end 
delay, and energy consumption. 

The multiobjective routing optimization problem described 
in the following can be solved using existing multiobjective 
search techniques |9|. However, the description of such a 
heuristic is out of the scope of this paper and will be addressed 
in later works. The paper is organized as follows. In Section [ll| 
we present our cross-layer framework based on a probabilistic 



network model. Section III formulates routing in an ad hoc 
network as a multiobjective optimization problem and Section 



Node Probability 

o 



Link Probability 



Node / Py fphys^ij) Node j 

Fig. \. Node and link probabilities on a link («, j). 
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TV| provides a first formulation applied to sensor networks. 
Results for a simple problem instance are then given in 
Section |V] to illustrate our modeling framework and Section 
VI concludes the paper. 



II. A Cross-layer Framework for Network 
Modeling 

A. Probabilistic network model 

Our proposed model considers a probabilistic network 
which is characterized by two probability measures: link and 
node probability. These two parameters completely character- 
ize the network and capture cross-layer interactions. 

The node probability {xd captures the availability of 
node i for routing purposes, i.e. the probability that node 
i re-broadcasts a received packet. The node probability has 
two components (xi = ■ Xi), one that is determined by 
the environment and protocol implementations at adjacent 
layers, ^i, (e.g. congestion models, node failures, security 
risks, energy levels), and one component Xi that coiTesponds 
to network routing choices, which we aim to optimize in the 
multiobjective routing framework. 

The link probability (p^) captures the link availability, 
i.e., the probability of a successful transmission over a link 
Characterization of the link probability is impacted by 
impairments and enhancements at various layers of the proto- 
col stack such as fading at the physical layer or congestion at 
the MAC layer. Both node and link probabilities are illustrated 
in Fig. [1] 

Both node and link probability measures are strongly related 
due to the nature of the wireless channel. Hence, once the 
node probabilities Xi ^re set, the activity of every node of 
the network is fixed and the interference distribution can 
be completely determined given the nodes 's activity on the 
wireless channel. As a consequence, the link probabilities 
can be computed as a function of the signal to noise and 
interference ratio (SINR). Once link and node probabilities 
are available, various performance metrics such as delay, ro- 
bustness or energy consumption can be calculated for various 
transmission schemes (unicast, multicast, broadcast, anycast, 
etc.). 

In the following, we consider the set of node probabilities as 
the variables of the network optimization problem. Finding the 
best possible routing with respect to one particular criterion 
reduces to the problem of selecting the set of node probabilities 
that optimizes one particular objective of the network. Within 
a multiobjective perspective, solving the network optimization 
problem requires finding the set of Pareto-optimal solutions 



that concurrently optimizes several performance metrics of the 
network. 

To illustrate our framework, we consider here a network 
where the nodes are independent and randomly distributed 
according to a random point process of density p over a disk 
T). The communication between any two nodes is performed 
in a half-duplex mode over a single to multi-hop path. The 
bandwidth of the channel is divided into R resources (time 
slots, frequencies or codes). For clarity purposes, we present 
this model in the context of time-multiplexing. 

This paper concentrates on a single flow but our framework 
can be extended to multiple flows since the proposed interfer- 
ence model accurately accounts for all the nodes transmitting 
in the network. Hence, one source transmits a constant traffic 
in one of the R time slots. A relay does not keep track of the 
packets already transmitted and consequently may forward the 
same packet several times. However, a node relays the packets 
in the order they are received in one of its available resources. 
If several packets are received in the same frame it can 
only transmit the proportion of packets its global transmission 
probability Xi allows. The packets that the node can not 
forward are dropped. The maximum number of hops Hm a 
packet can travel in the network is also fixed. 

B. Link probabilities 

A realistic link in time slot r is characterized by 

its transmission probability Pij{r), which is a function of 
the statistical distribution of the SINR at the location of the 
destination node j. Such a computation captures the cross- 
layer impact of the routing decision on the physical layer 
performance since the activity of all the nodes of the network 
are accounted for statistically in the model. The following are 
some preliminary definitions and notations that are needed to 
define the link probability: 



Pathloss attenuation factor: 



reflects the attenuation 



due to propagation effects between node i and j. In our simu- 
lations, the simple isotropic propagation model is considered. 

Interference: Since we consider time-multiplexed chan- 
nels, interference only occurs between transmissions using 
the same channel at the same time. Hence, the power of 
interference Iij{r) on a link using resource r and 

computed at node j is defined by: 



K 



(1) 



k=l 



where K is the number of interfering signals in resource r. 

SINR: The SINR between any two nodes i and j in 
resource r is given by: 



(2) 



No + hj{r) 



where Pij is the power received in j, Iij (r) is the interference 
power on the link and A^g the noise power density. We have 
Pij = Pi 0-ij for a fixed nominal transmission power Pi and 
a pathloss attenuation factor a^ . 



Packet error rate (PER): For a specific value of SINR 
7, the packet error rate PER can be computed according to: 

PER{-f) = 1 - [1 - BER{-/)f'' (3) 

where Nf, is the number of bits of a data packet and BER{"i) 
is the bit error rate for the specified SINR per bit 7 which 
depends on the physical layer technology and the statistics 
of the channel. Results are given for an AWGN channel 
and a BPSK modulation without coding where BER{'^) ~ 
Q (^/2^) = 0.5 * erfc(/7). 

Transmission rate: The activity of a network node in a 
channel re [1, .., i?] is given by its transmission rate Ti{r) E 
[0, 1] in that particular channel. This rate is defined as the 
percentage of time a node i transmits using resource r. 

Additional Notations: A node i is said to be active in 
the network if '^i ('') > 0, and 

- M gives the number of active nodes of the network, 

- An interfering set on a link (i, j) is a set of X < M — 1 
active nodes, 

- C-i refers to the set of all possible interfering sets and 
has a cardinality of i = T.k=i (^\^^) + 1- 

The link probability: Pij{r) depends on the distribution 
of the SINR, and consequently on the distribution of the 
corresponding packet error rates. It is defined by the equation: 

L 
1 = 1 

where the index I represents one of the L interfering sets. 
Consequently, 7;(r) is the SINR experienced because of 
the interfering set I on the link for the resource r 

and PERi{r) is the corresponding PER. The SINR can be 
computed according to Eq. (|2| considering the K interfering 
links of / and the PER according to Eq. ([3]). 

P;(r) is the probability for the link to experience the 
interference distribution I in resource r, i.e. the probability that 
the nodes of the interfering set I are transmitting concurrently 
and the others are not. Hence, this probability for a link 
is given by: 

K M-K-l 
k=l 7n=l 

In Eq. (|5]), Y[k=i'''k(T) gives the probability that the K 
active nodes of the interfering set I are transmitting and 
nm=f ""^1 - Trn{r)) the probabiHty that the M - K - I 
other active nodes are not. 

C. Node probabilities and transmission rate 

The variables of our model are the probability Xi ~ for 
each node i to re-transmit a received message. In the following, 
we consider that = 1 to simplify our model. Hence, the 
main variable is the 'forwarding probability' Xi. There is no 
notion of routing paths herein and a packet sent by a source 
may use one or more paths in parallel to reach the destination. 
For Xi = 1 each received packet by node i is forwarded. For 



Xi <1 node i drops the packets with probability 1 — Xi. Values 
of Xi E ]1, i?] are not allowed yet as they imply that node i 
transmits several copies of the same packet. 

As stated earlier, the transmission rate T,;(r) in resource 
r is a function of the node probability Xi but also depends 
on the amount of traffic coming into node i, which is a 
function of the activity of the other nodes of the network. 
As a consequence, computing the values of (r) knowing the 
Xi values is intractable since determining the Ti{r) requires 
the knowledge of the link probabilities which are themselves 
a function of the Ti{r) values. However, the reverse approach 
where the variables x are expressed as a function of the Ti{r) 
can be easily derived as stated below. Hence, such a reverse 
approach leads to the use of the transmission rates as the 
variables of our multiobjective optimization problem instead of 
the forwarding probabilities. This reverse approach represents 
an important contribution of our cross-layer model since it 
captures an exact picture of the interference distribution at 
the physical layer and determines the corresponding node 
forwarding probability Xi at the routing level. 

Relationship between Xi and the Ti{r): Given the 
values of Ti{r),yr e e [l..iV], we can define the 

quantity of information coming from all the neighbors of node 
i (except from the destination) by: 

Qi^ ^ ^Pkt{r).Tk{r).Vki (6) 

k^{i,D} r 

where Pki{i')-Tk{f)-Vki is the probability that a packet arrives 
in node i from node k in resource r. 

The variable Vki is introduced to represent the usefulness 
of the link (k,i) with respect to the maximum number of 
hops constraint. Hence, if no data can arrive from neighbor 
k because the hop count h for all the packets k received is 
already equal to Hm, we have Vki = 0. On the contrary, we 
have Vki — \ if k only receives packets with a number of 
hops h < Hm- If k receives packets with both h < Hm 
and h — Hm, v^i represents the proportion of packets being 
retransmitted. 

The quantity of information going out of i is given by 
the sum of the Ti{r) over all the time slots. Hence, we can 
determine the global forwarding probability of i to be: 

= ^r^jr) 

J2k^{i,D}J2rPki{r)■Tk{r).Vk^ 
III. A MULTIOBJECTIVE OPTIMIZATION PROBLEM 

The performance of most wireless networks can be assessed 
with regards to various criteria such as throughput or capacity, 
end-to-end transmission delay, overall energy consumption or 
transmission robustness. The purpose of the multiobjective 
framework presented in this work is to determine, given a 
network and a communication pattern, what kind of trade- 
offs arise between chosen performance metrics when varying 
the routing strategies. It relies on the cross-layer probabilistic 
network model presented in Section |ll] 



A. Variables of the Multiobjective (MO) Framework 

The routing strategies are the variables of our muhiobjective 
optimization problem and a solution is defined by: 

Definition 1 A solution S of the MO framework is defined 
by the set of transmission rates Tj(r) G [0,1] used by each 
node i on each resource r: 



S^{n{r)} 



ie[l,...N],re[l..R] 



(8) 



The set of node probabilities jy] is derived according 

to Eq.(|7| and represents the routing strategy of the network. 
Each variable T,;(r) takes its values in a discrete set F of size 
T = |r|. As a consequence, the solution space is derived as: 
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( V) T 



R.m 



(9) 



In order to reduce the size of this very big search space, 
we only consider solutions where at least one cumulative 
time slot per node is available in the frame, i.e. s.t. Vi G 
[l,iV], X^tLi '^ii''') < ^ — 1- The solutions that do not meet 
this constraint are usually very bad solutions since at least one 
of the nodes of the solution is transmitting in all its time slots 
preventing a failure free packet reception. 

Using this definition of a routing strategy, a solution may 
reflect various features: it can be single-hop or multi-hop, 
single path or multi-path, probabilistic or deterministic. The 
aim of our MO framework is to obtain the set of Pareto-optimal 
routing strategies of the MO problem. A Pareto-optimal set 
is composed of all the non-dominated solutions of the MO 
problem with respect to the performance metrics considered. 
A solution A dominates a solution B for a rt— objective MO 
problem if A is at least as good as B for all the objectives 
and A is strictly better than B for at least one objective. 

IV. A FIRST APPLICATION TO SENSOR NETWORKS 

We propose in the following to assess the performance 
of a wireless sensor network (WSN) by capturing the trade- 
offs that arise between end-to-end robustness, overall energy 
consumption and end-to-end delay. These criteria are prevalent 
since providing a maximal network throughput is usually not 
the main task of a WSN. The criteria are defined for a single 
source-destination pair (5, D). 

A. Robustness criterion 

Robustness is defined as the probability that a message 
emitted at S successfully arrives at D in at most Hm hops. 
The robustness criterion is given by: 



(10) 



For any two nodes i and j of the network, T"^ represents the 
event that a message transmitted by i successfully arrives in 
j in at most H hops. Our aim is to maximize V{Tg^'). 
Definition 2: Global Unk probability. 



For a link the global link probability pij is the 

probability that a message arrives with success at node j. It 
is given by: 



R 

E 



P^l= 2^ P^Ar) v"'r-rr) ^^^^ 

where Pij{r) is the link probability between i and j for 
resource r (cf. Eq. Q) and Ti{r) / J^r '''ii''') the probability 
for the packet to be sent using r. 

Definition 3: Robustness probability. 

'P[Tg^') is the probability that the message arrives success- 
fully in D in at most Hm hops and is given by: 

'P{T^^') = l-Y[{^-'P{TsD\H = h)) (12) 

where V{Tsd\H = h) is the probability for a packet to 
arrive in h hops at D. For h = 1, 'P{Tsd\H = 1) = Psd, 
the successful transmission probability on the link {S, D) 
For /i > 1, we have: 



following Eq. ( 1 1 

V{TsD\H^h) 



1 



Ns 

nil 



Psj Xj V{Tjd\H 



-1)] 
(13) 

with Ns the number of possible first hop relays of S\ psj the 
link probability between S and its neighbor j; V{TjD\H = 
h — 1) the probability to reach D in [h — 1) hops and Xj the 
forwarding probability of j. The set of Ns relays is given by 
all the nodes different from S that are active in at least one of 
the time slots in the current solution (i.e. having X]tLi(^i) > 
0, ^ ^ {j, S)). 

To reduce the computation complexity of the robustness 
probability, a restricted set Ns of first hop relays may be 
considered but the loss in terms of accuracy is hard to quantify. 
Therefore, we rather introduce a link threshold value Vth 
computed for each path made of h hops. While recursively 
calculating V{Tsd\H = h), if the probability of a path gets 
lower than Vth^ the recursion is stopped for that particular 
path and its contribution to V{Tsd\H = h) is set to zero. 

B. Delay criterion 

The end-to-end delay is the sum of the times spent at each 
relay on a multi-hop path where each relay introduces a delay 
of 1. The criterion fo is defined by: 



fD = R- 



j2ih~ir.R, 



(14) 



The quantity {h— 1) is the delay needed by a packet to arrive 
in h hops using {h — 1) relay nodes. The scaling factor R 
represents the delay induced by the R resources. Rk is the 
probability that the packet arrived in exactly h hops and did 
not arrive in 1, or 2 or (/i — 1) hops. For h = 1, we have 
R^ = P{TsD\h ^ 1) and for h > I: 

h-l 

Rh = V{Tsd\H = h). [](1 - V{Tsd\H = i)) (15) 



Transmission Power 151mW Nq -154dBm/Hz 

Bandwidtli 1Mbps / 2AGHz 

Pathloss exponent a 3 Cliannel Model AWGN 

Antenna gains Gt=Gji=1 Modulation BPSK 



Fig. 2. Propagation and physical layer parameter values. 

If no route exists between S and D then fo = +oo. 

C. Energy criterion 

The energy criterion Je is given by the total forwarding 
energy needed for a packet sent by S to reach D. We do not 
account for the energy spent by the initial transmission in S. 
The reception (resp. transmission) of a packet at node j in 
resource r consumes ef{r) (resp. ej{r)). Hence, the energy 
criterion is defined as: 

fE^J2^(^SD\H = h) (16) 

h=l 

where 8{Tsd\H — h) is the total energy needed by the /i-hop 
communications between S and D defined by: 

£{TsD\H = h) = 

Ns 

J2 {Psyef+psyx,. [ej + £{T,d\H ^h-l)]) (17) 

In Eq. ([17]), psj.ef is the energy consumed for a packet 
reception by the neighbor j of S; psj-Xj.eJ is the energy 
consumed for the packet transmitted by neighbor j and 
psj-Xj.£{TjD\H = h — 1) is the total energy consumed by 
the following possible paths made of (/i — 1) hops between 
neighbor j and the destination. For h = 1, £{Tsd\H = 1) = 
since the energy in S is not accounted for 

V. First Results 
A. M-Relay problem 

The results presented in this section are obtained for a 
small problem instance for two reasons. First, we are able 
to determine the complete Pareto-optimal set of solutions 
using an exhaustive search. Secondly, such a problem can 
be easily analyzed and provides a first illustration of our 
multiobjective framework. We will tackle bigger instances 
using multiobjective optimization algorithms |9|. 

In the following, the network is composed of iV = 333 
nodes uniformly distributed with density p — 0.004 over a disk 
T) of radius The distance between S and D is of about 
215 meters. To reduce border effects, S and D are selected 
within a radius Rq « R-d which ensures that the power of 
a node at distance Rc is below the noise power for the nodes 
located at distance Rj). We consider R ~ 2 time slots and use 
a probabilistic discrete variable space where Ti{r) takes its 
values in the set T = {0, 0.05, 0.1, . . . 0.9, 0.95, 1.0} of |r| = 
21 elements. A link robustness threshold of Vth — 10"^° is 
set. Propagation and physical layer parameters are summarized 
in Fig. [2] 



The dimension of the search space can be modified by 
setting a maximum number of forwarding nodes M in a 
solution S. This sub-problem is addressed in the following 
as the Af -relay problem instance. 

B. Pareto-optimal set for the 1 -relay problem 

In this problem instance, we set AI = 1 and Hm = 2. 
In that particular case, the search space has a dimension of 
72813 solutions and the Pareto-optimal set is obtained with 
an exhaustive search. 

For this instance, the direct link {S, D) is very weak. A 
robustness of only V{Tg^') — 0.0003 is achieved with a delay 
of = and an energy of fs = 0. Only 24820 solutions 
fulfill the constraint Xi < 1 that forbids a node to duplicate 
packets. Among these solutions, 3855 solutions are Pareto- 
optimal, representing respectively about 5% and 15% of the 
whole and the constrained solution space. For all the Pareto- 
optimal solutions the relay never transmits in the first time slot 
concurrently with the source. Hence, the model suggests that 
the Pareto-optimal set for this case is composed of solutions 
that minimize interference. 

The performance of the Pareto-optimal set of solutions 
is represented in Fig. |3] in the space defined by the three 
evaluation functions. For clarity purposes, the projections of 
the Pareto set on the robustness/delay, robustness/energy and 
the delay/energy planes are displayed. The plots of Fig.[3]show 
that an improved robustness is obtained at the price of an 
increase in delay and energy. The trade-off between robustness 
and delay can be easily understood since higher robustness is 
achieved when the relay contributes with a higher forwarding 
probability Xi, inducing an increase in delay. Similarly, an 
increase of Xi triggers an accrued average energy consumption 
since the relay is forwarding packets more often. 

The energy consumption for all the Pareto-optimal solutions 
belongs to a discrete set of 21 energy levels which is a direct 
consequence of the 21 values of Ti{r) defined in this problem 
instance. Hence, the definition of a continuous transmission 
rate variable Ti{r) would provide the most precise description 
of the Pareto set. However, tackling the continuous formulation 
of our problem is much more challenging and for our study, 
we will stick to the simpler discrete formulation which still 
provides a fair representation of the Pareto set. 

The Pareto set is composed of solutions where relays belong 
to a set of 226 nodes, which represents about two thirds of the 
number of nodes of the network. The location in the network of 
these 226 nodes is presented in Fig. [4] We also highlighted on 
this figure the relays that provide a near perfect transmission. 
We can conclude that the relays located in an ellipse near the 
middle of the (S*, D) distance provide the best robustness at 
the price of the highest delay and energy. The other relays 
present in the Pareto set provide various trade-offs depending 
on their values of T,;(r). 

This first simple study shows that the proposed multiob- 
jective probabilistic network model provides a coherent and 
complete view of the trade-offs that arise between robustness, 
delay, and energy in our network. A more extensive analysis 



Fig. 3. Representation of the projections of the Pareto-optimal set for the 1 -relay problem. 
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model. Future work will extend the model to consider more 
complex network scenarios, such as to account for various 
network topologies, to consider multiple concurrent flows in 
the network, and to use more refined cross-layer interactions 
and interference models. 
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Fig. 4. Location of the nodes that provide Pareto-optimal solutions in the 
network (blue crosses) and of the nodes that provide a near quasi-perfect 
robustness (full black dots), i.e. fn > 0.999. 



of the performance of the model has to be performed next 
by considering a complete solution space and various network 
topologies. For such instances, the size of the search space 
prohibits the use of an exhaustive search of the Pareto set. 
Hence, we will concentrate on implementing a combinatorial 
multiobjective optimization algorithm to obtain the best pos- 
sible representation of the Pareto-optimal set. 

VI. Conclusion 

In this paper, we have proposed a novel multiobjective 
optimization framework for network routing in wireless ad 
hoc networks. Our proposed framework consists of a gen- 
eral probabilistic network model capable of capturing the 
impact and interaction of a wide range of resource/interference 
management techniques, and various channel conditions and 
network scenarios. Used in conjunction with metaheuristic 
optimization techniques, this framework provides an efficient 
tool to capture the trade-offs between different performance 
metrics and obtain bounds on the achievable performance 
of routing for a single source-destination transmission. Pre- 
liminary results were obtained in characterizing the delay, 
robustness, and energy tradeoffs for a 2-hop sensor network 
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